Tensor-force effects on single-particle levels and proton bubble structure around the Z 

or N = 20 magic number 
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Applying the semi-realistic AW interactions that include realistic tensor force to the Hartree-Fock 
calculations, we investigate tensor-force effects on the single-particle levels in the Ca isotopes, ft 
is clarified that the tensor force is important in reproducing the experimental difference between 
e(plsi/ 2 ) and £(p0d 3 / 2 ) (denoted by Aeu) both at 40 Ca and 48 Ca. The tensor force plays a role 
in the A^-dependence of Aeis also in neutron-rich Ca nuclei. We further investigate possibility 
of proton bubble structure in Ar, which is suggested by the pls 1 / 2 -p0d 3/ / 2 inversion in 4S Ca and 
more neutron-rich Ca nuclei, by the spherical Hartree-Fock-Bogolyubov calculations. Even with the 
inversion at 48 Ca the pair correlation prohibits prominent bubble distribution in 46 Ar. Bubble in Ar 
is unlikely also near neutron drip line because either of unboundness or of deformation. However, 
34 Si remains a candidate for proton bubble structure, owing to large shell gap between pls x / 2 and 
p0d 5/2 . 

PACS numbers: 21.10.Pc, 21.10.Ft, 21.60.Jz, 27.40.-fz, 27.30. +t 
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Introduction. Owing to the progress in experiments on 
unstable nuclei, it has been recognized [l| that the nu- 
clear shell structure depends on Z or N as often called 
shell evolution. Moreover, it is now known that the 
tensor force, which should be contained in the nucleon- 
nucleon (NN) interaction, plays a crucial role in the shell 
evolution. While the tensor force had been ignored in 
the conventional mean-field (MF) or the energy-density- 
functional (EDF) approaches, there have been several at- 
tempts to incorporate the tensor force into the those ap- 
proaches; e.g. the calculations with the Skyrme [5|-|8( 
or the Gogny [9J] interactions. However, without well- 
established strengths and/or radial-dependence of the 
tensor force, it is not straightforward to pin down tensor- 
force effects quantitatively from those interactions. One 
of the authors (H.N.) has developed the so-called semi- 
realistic NN interactions [ltj, which is applicable to the 
MF calculations. The recent parameter-sets, M3Y-Pn 
with n > 5 dH [ll| , include the tensor force originat- 
ing from the G-matrix at the nuclear surface. Because 
of the realistic nature of the tensor force in them, these 
semi-realistic interactions are suitable to investigate the 
tensor-force effects on the shell evolution. 

Although the MF or EDF approaches give single- 
particle (s.p.) levels in a self-consistent manner, each of 
the s.p. levels does not correspond to a single observed 
state even beside the doubly magic nuclei, being frag- 
mented over a certain energy range. The s.p. energies in 
the MF approaches should better correspond to the av- 
eraged energy of the states having the specific quantum 
numbers weighted by the spectroscopic factors. There are 
not many cases in which the spectroscopic factors have 



exhaustively been measured with good accuracy. The 
proton lsi/2 and 0d 3 / 2 hole states in 40 Ca and 48 Ca pro- 
vide us with indispensable examples, for which sum of 
the measured spectroscopic factor exceeds 90% 12. 13j|. 
These experimental data have disclosed a notable con- 
sequence that the plsj/2 and p0d 3 / 2 levels invert as N 
increases from 40 Ca to 48 Ca. The plsi/ 2 -p0d 3 / 2 differ- 
ence and their inversion seem to supply a good test of the 
MF effective interactions (or EDFs). Since it has been 
suggested that the tensor force plays a important role in 
this inversion 1J, [l5j , it is of interest to apply the semi- 
realistic interaction. Moreover, the pls 1 / 2 -p0d 3 / 2 inver- 
sion at 48 Ca suggests a possibility of proton "bubble" 
structure, depletion of the proton density at the center 
of the nucleus, in the two-more-proton deficient nucleus 
46 Ar It is noted that similar inversion has been 

predicted in several Ca nuclei near the neutron drip line, 
which suggests proton bubble in Ar. Possibility of proton 
bubble structure has been pointed out also for 34 Si [l7j |. 
Although it has been difficult to measure charge densities 
at the center of unstable nuclei, the new technology such 
as SCRIT 18j could open a way to observe such proton 
bubble structure. 

In this paper we apply the self-consistent Hartree-Fock 
(HF) and Hartree-Fock-Bogolyubov (HFB) calculations 
with the semi-realistic interactions to the Si to Ca nuclei, 
and investigate tensor-force effects on the shell evolution 
and the bubble structure. 

Effective Hamiltonian. We apply the spherical HF 
and HFB calculations by using the Gaussian expan- 
sion method The details of the basis func- 
tions are given in Ref. [HI- The effective Hamiltonian 
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has the form H = Hjv + Vc — H c .m., where Hn{= 
2,pf/2M + ^Z^jVij), V c and H c , m , denote the effec- 
tive nuclear Hamiltonian, the Coulomb interaction and 
the center-of-mass Hamiltonian, respectively. It is noted 
that the exchange term of Vc is treated exactly and that 
both the one- and the two-body terms of H c . m , are sub- 
tracted before iteration. 

The M3Y-Pn semi-realistic interactions have been ob- 
tained by modifying the so-called M3Y-Paris interac- 
tion 22J, which was derived by fitting the Yukawa func- 
tions to the G-matrix at the nuclear surface. Density- 
dependent contact terms have been added to the M3Y- 
Paris interaction so as to realize the saturation, and the 
LS channels have been enhanced in order to reproduce 
the is splitting at the MF level. We employ the M3Y- 
P5' [23J and M3Y-P7 11| parameter-sets in the following. 
In both sets the tensor force of the M3Y-Paris interac- 
tion is maintained without any change, as in M3Y-P5 [3 1 
and P6 11 1. For comparison we also use the DIM [24 1 
parameter-set of the Gogny interaction, which has no ten- 
sor channels. Whereas we have implemented calculations 
with M3Y-P6 and D1S [H, results of M3Y-P6 (D1S) are 
similar to those of M3Y-P7 (DIM). 

TV -dependence of plsi/2~p0d 3 /2 levels in Z = 20 nu- 
clei. We here express the s.p. energy difference under 
interest as Aei 3 = e(plsi/ 2 ) — e(p0d 3 / 2 ). Figure Q] shows 
TV-dependence of Aei3 in the Ca isotopes obtained by 
the spherical HF calculations. The experimental values 
of Ae\ 3 in 40 Ca and 48 Ca are also displayed for compar- 
ison, which are obtained after average weighted by the 
spectroscopic factors [12, (HI- Although significant TV- 
dependence is found in the experiments on Ae\ 3 of the 
Ca isotopes, not many MF interactions (or EDFs) can re- 
produce this TV-dependence quantitatively [3]. In prac- 
tice, despite agreement of Ae\ 3 with the data at 48 Ca, 
there is significant discrepancy at 40 Ca in the DIM re- 
sult, as viewed in Fig. [TJ With D1S, Ae\ 3 is slightly 
shifted downward and the inversion at 48 Ca does not oc- 
cur. The slope of Ae 13 from 40 Ca to 48 Ca in the DIM 
result is typical to the MF calculations with no tensor 
force. On the contrary, as depicted in Fig. [TJ the semi- 
realistic M3Y-P5' and P7 interactions successfully repro- 
duce the TV-dependence of Ae 13 . In Fig. Q] we also show 
Z\e 13 in which contribution of the tensor force is removed 
from the M3Y-P7 result. Then Ae 13 varies in parallel to 
those of DIM, unable to reproduce the observed slope. 
This confirms the crucial role of the tensor force in Ae\ 3 
as nQfi/2 is occupied. 

The variation of Ae\ 3 from TV = 20 to 28 is a result 
of occupation of the n0f 7 / 2 orbit. As nO/7/2 is occupied, 
the tensor force tends to lower p0d 3 / 2 but not plsi/2, 
therefore increasing Aei 3 . In DIM the parameters in the 
central and LS channels have globally been adjusted to 
experimental data, and the inversion of the s.p. levels at 
48 Ca is reproduced. However, because lacking the ten- 
sor force, this inversion takes place at the expense of the 
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FIG. 1. (Color online) Asiz of the Ca isotopes. Green, blue 
and red lines represent the results with the DIM, M3Y-P5' 
and P7 interactions, respectively. Thin red dashed line is 
obtained from M3Y-P7 but by removing contribution of the 
tensor force. 



discrepancy in Z\ei 3 at 40 Ca. It should be commented 
here that there might be influence of ground-state cor- 
relations on Asi 3 which are not well taken into account 
in the averaged s.p. energies, and that its assessment is 
desired for complete understanding. 

The tensor force affects Aei 3 from TV = 34 to 40, i.e. 
as nO/5/2 is occupied. Experiments on the plsi/2 and 
the p0d 3 / 2 hole states around 60 Ca may provide further 
evidence of the tensor- force effects, if carried out in the 
future. Although it depends on the effective interactions 
even in MF approaches where the neutron drip line of Ca 
is, plsi/2~p0d 3 /2 inversion could occur again toward 70 Ca 
as n0<79/2 is occupied. However, whereas the inversion 
takes place already at 60 Ca with DIM, the inversion is 
delayed in the M3Y-P5' and P7 results, in which the 
parameters have been determined in the presence of the 
tensor force. 

Investigation of proton bubbles. The plsi/2-p0d 3 / 2 in- 
version in 48 Ca indicates that dominant configuration 
of the two-proton-deficient nucleus 46 Ar is (plsm) -2 - 
Since only the s-states give sizable density at the center 
of nuclei, 46 Ar is a candidate of a nucleus having proton 
bubble structure [l6j]. While it is not easy to measure 
matter or neutron densities in a mo del- independent man- 
ner particularly for unstable nuclei, charge densities may 
unambiguously be extracted from the electron scattering 
experiments. Since the M3Y-Pn semi-realistic interac- 
tions successfully describe the TV-dependence of Ae\ 3l in- 
dicating that mechanism giving rise to the plsi/2-p0d 3 /2 
inversion is correctly contained, it will be of certain inter- 
est to investigate possibility of proton bubble structure 
by applying these interactions. 

We present proton density distributions at 48 Ca and 
46 Ar obtained from the MF calculations with DIM and 
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and accordingly the proton bubble structure, as indicated 
by Fig. [TJ As previously discussed this is mostly due 
to the occupation of nO/5/2. In contrast, the inversion 
is predicted by all the interactions under investigation 
around 70 Ca. The 70 Ca nucleus is bound in the HFB 
calculations with DIM, M3Y-P5' and P7 
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FIG. 2. (Color online) Proton density distributions in 48 Ca 
and 46 Ar obtained from HF and HFB calculations. 



M3Y-P7 in Fig. O The pair correlation is quenched 
in the ground-state of 48 Ca and the HF and HFB re- 
sults are identical. Within the spherical HF regime with 
M3Y-P7, we have depletion of the proton distribution 
around the origin at Ar, since plsi/ 2 becomes unoccu- 
pied if compared with 48 Ca. The same holds for M3Y- 
P5', though not displayed. Such depletion is not found in 
the HF result with DIM, in which the ground state has 
the (j>0d 3 / 2 )~ 2 configuration. Despite the plsi/z-pOda/z 
inversion at 48 Ca, the total energy of the (p0d 3 / 2 )~ 2 
state becomes lower than that of the (plsi/ 2 ) -2 state 
at 46 Ar in the DIM result. On the other hand, once 
the pair correlation is taken into account, small energy 
difference between pls\/ 2 and p0d 3 / 2 significantly mixes 
up the (plsi/ 2 )~ 2 and the (p0d 3 / 2 )~ 2 configurations for 
both interactions. Thus the depletion observed in the HF 
densities with M3Y-P7 is smoothed out when the pairing 
is switched on. 

While one might consider the formation of proton bub- 
bles around 60 Ca because of the pls 1 / 2 -p0d 3 / 2 inversion 
given by DIM, the tensor force prohibits the inversion 
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a glance, the plsx/ 2 -pQd 3 / 2 inversion exhibited in Fig. Q] 
seems to suggest proton bubble structure in 64 ~ 68 Ar even 
with M3Y-P7. However, the 64 - 68 Ar nuclei lie beyond 
the neutron-drip line while 66_70 Ca are bound, within 
the spherical HFB calculation. Although these Ar nu- 
clei might be deformed and thereby their energy could 
be lowered enough to be bound, the deformation signifi- 
cantly mixes the (plsi/ 2 ) -2 state with others, which eas- 
ily destroys bubble structure. We therefore conclude that 
proton bubble structure is unlikely to be observed in any 
of the Ar nuclei. 

While the N = 20 magic number disappears in the 
proton-deficient region of Z < 12, it remains magic in 
Z > 14, keeping the nuclear shape spherical. Possibil- 
ity of proton bubble structure has been pointed out also 
for 34 Si The energy difference Aei$ = e(pls 1 / 2 ) — 
e(p0d 5 / 2 ) exceeds 5MeV for 36 S in the spherical HF cal- 
culations, irrespective of the effective interactions under 
consideration. We note that — zAei 3 is less than 2MeV, 
giving rise to sizable pair excitation at 36 S. The occupa- 
tion probability of plsi/ 2 is 0.6 — 0.7 in the spherical HFB 
results. However, owing to the large subshell gap Zieis, 
the ground state of 34 Si is expected to be (plsi/ 2 ) -2 with 
good approximation. Proton density distribution of 34 Si 
is depicted in Fig.[3]in comparison with that of 36 S. Since 
Z\ei5 is sufficiently large, prominent proton bubble struc- 
ture is predicted in 34 Si in the MF calculations for all 
the effective interactions considered in this paper. The 
large Aei§ prevents the pair correlation from mixing the 
{plsi/ 2 )~ 2 state with the other states {e.g. (p0g? 5 / 2 )~ 2 ) 
in the ground state of 34 Si, giving identical density distri- 
bution between HF and HFB. Although influence beyond 
the MF regime is not obvious [jjj ]. 34 Si remains to be a 
candidate for the proton bubble. Future experiments on 
the charge density of this nucleus are awaited. 

Summary. We have investigated tensor-force effects on 
the shell evolution and the bubble structure around the Z 
or N — 20 magic number, by applying the self-consistent 
HF and HFB calculations with the semi-realistic inter- 
actions. As well as having confirmed significance of the 
tensor force in the A-dependence of the plsi/ 2 -p0d 3 / 2 
difference in the Ca isotopes, we have shown that realis- 
tic tensor force gives adequate A-dependence if combined 
with appropriate central and LS forces, as in the M3Y-Pn 
semi-realistic interactions. Although the plsi/ 2 -p0d 3 / 2 
inversion is predicted for the Ca nuclei near the neutron 
drip line, the inversion is delayed with the interactions 
including realistic tensor force, suggesting the inversion 
only in 66 - 70 Ca. 

On the basis that semi-realistic interactions reproduce 
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FIG. 3. (Color online) Proton density distributions in 36 S and 
34 Si obtained from HF and HFB calculations. 



the observed plsi^-pOc^^ inversion reasonably well, we 
apply them to perform predictions for the proton-bubblc- 
structure problem. In Ar, although we view depletion 
of the proton density within HF, it is predicted that pair 
correlation prohibits the bubble structure from being re- 
alized. Although one may anticipate proton bubble struc- 
ture in the highly neutron-excess nuclei 64-68 Ar because 
of the inversion in 66 ~ 70 Ca, these nuclei might be un- 
bound or deformed, and therefore proton bubble struc- 
ture is unlikely. On the contrary, it may remain hopeful 
that the proton bubble structure will be observed in fu- 
ture experiments in 3 Si. 
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